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The lightest Kaluza-Klein particle appearing in models with universal extra dimensions has re- 
cently been proposed as a viable dark matter candidate when the extra dimensions are compactified 
on a scale of the order of 1 TeV. Underlying assumptions of this proposal are that the size of 
the extra dimensions stays constant and that the evolution of the universe is given by standard 
cosmology. Here we investigate, both analytically and numerically, whether this is possible with- 
out introducing an explicit stabilization mechanism. By analyzing Einstein's field equations for 
a (3+n + 1) - dimensional homogeneous, but in general anisotropic universe, we find that approxi- 
mately static extra dimensions arise naturally during radiation domination. For matter domination, 
however, there are no solutions to the field equations that allow static extra dimensions or the usual 
behavior of the scale factor for ordinary three-dimensional space. We conclude that an explicit 
mechanism is needed in order to stabilize the extra dimensions and reproduce standard cosmology 
as we know it. 



I. INTRODUCTION 



Everyday experience seems to suggest that our world consists of four space-time dimensions. However, at the 
beginning of the 20 th century Nordstrom, Kaluza, and Klein (KK) already realized that this may in fact not be the 
case In the last few years there has again been a great deal of interest in models with extra dimensions, most 

notably due to the influence of string (or M) theory which in its usual formulation requires more than four dimensions 
(see e.g. Q). In particular, this has led to a number of brane- world scenarios where all, or only some, of the gauge 
bosons are allowed to propagate in the extra dimensions while matter fields are restricted to (3 + 1) - dimensional 
branes. For a review of different models with extra dimensions, see, for example, || and references therein. 

In this context, a specific model of so-called universal extra dimensions (UED) has recently been proposed by 
Appelquist et al. [||, in which all standard model fields are allowed to propagate in the extra dimensions. As usual, 
quantization of the extra-dimensional momentum leads to a tower of KK states that appear as new massive particles 
in the effective four-dimensional theory. The existing constraints on electroweak observables translate into bounds on 
the compactification scale R, which is related to the mass of the lowest excitations by M ~ 1/R. For one or two UED 
these bounds are of the order of a few hundred GeV and thus within reach of the Large Hadron Collider (LHC) or 
the Fermilab Tevatron run II [E 0, H 13 . 

The UED model is not only of great interest from the point of view of particle physics [H HJ 13 El, 

but might 

also provide a solution to one of the most outstanding puzzles in modern cosmology - the nature of dark matter 
[Til flH H3|. In the UED scenario the lightest KK particle (LKP) is stable because of KK parity conservation and 
could therefore still be present today as a thermal relic. Furthermore, if it is also neutral and nonbaryonic it has all 
the properties of a weakly interacting massive particle (WIMP), one of the most promising dark matter candidates 
(see E3 f° r a nice introduction to WIMP dark matter). According to [Iq. both the KK photon (B^) and the KK 
neutrino could account for dark matter with Q,m ~ 0.3, as suggested by the current cosmological concordance model 
[l8|. if one assumes a compactification scale of about R ~ 1 TeV -1 size. Indirect and direct detection properties of 
such KK dark matter candidates are promising for the next generation detectors [l^. I20L l2ll |22|. In fact, the KK 
neutrino seems to be ruled out already by present data [2lJ . 

However, when considering the freeze-out process and the further cosmological evolution of thermally produced 
LKPs, the size of the extra dimensions has so far been assumed to stay constant - although no explicit stabilization 
mechanism has been given. As already noted in |l6l |23| , the resulting relic density today depends crucially on this 
assumption, and it is therefore important to investigate whether it can be justified within the UED framework. Since 
the UED may be relatively large (in fact, they must be if the LKPs are not to overdose the universe), their evolution 
should be governed by Einstein's field equations. The aim of this work is therefore to carefully study the dynamics 
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of the universe as described by an appropriate extension of the usual Friedmann equations that results from the 
field equations in higher dimensions. Specifically, we focus on solutions with constant or only slowly varying extra 
dimensions in the absence of any explicit stabilization mechanism. However, our analysis is more general in that we 
investigate whether there is any time evolution of the extra dimensions that corresponds to standard cosmology in 
four dimensions; in this context, we also include a numerical study of the transition regime between the radiation- and 
matter-dominated eras. For earlier work on the stability properties of hi ghe r-dimensional cosmolo gies , often referred 
to as Kaluza-Klein or multidimensional cosmology, see, for example, [24T.I25L I2H I27I I2H , 12^. I30L l3lL l32 fl and references 
therein. 

This paper is organized as follows: In Sec. [H] we introduce the cosmological solutions to Einstein's field equations 
for a (3 + n+ 1)- dimensional homogeneous, but in general anisotropic, universe. Here, we also comment on the 
interpretation of pressure in higher dimensions and derive a general relation between pressure and energy density in 
UED cosmology. Necessary conditions that every solution with static extra dimensions in such a model must satisfy 
are then derived in Sec. IIIII In the next two sections we study the existence of solutions with (nearly) constant extra 
dimensions during radiation and matter domination, respectively. A possible transition between these two regimes is 
then outlined in Sec. IVII Finally, Sec. IVIll discusses the implications for the UED scenario and concludes. 

II. SETUP 
A. Basic equations 

We introduce n universal extra dimensions and adopt coordinates X A , A = 0, 1, . . . , 3 + n, with 

x» = X» (p = 0,1, 2, 3) (1) 

and 

x l = X l (i = 1,2,3) (2) 
being the coordinates for ordinary four-dimensional spacetime and three-dimensional (3D) space respectively, and 



y 



p — 



X 3 +p (p=l,...,n) (3) 



the coordinates for the UED. In the absence of a cosmological constant, Einstein's field equations are given by 



G A B = R A B -\R8 A B = K*T A B . (4) 



Here, n 2 is defined as 



K 2 



M 2+n ' 



(5) 



where M is the higher-dimensional Planck mass. In the case of compactified extra dimensions with volume Vr n ) it is 
related to the usual Planck mass by |3^| 

M 2 Pl =V {n) M 2+n . (6) 

In the UED scenario, there is no localization mechanism that confines particles to a brane, so we assume that 
the energy density is distributed homogeneously throughout all dimensions. Thus, we are looking for homogeneous 
solutions to the field equations that are isotropic in ordinary three-dimensional space and also - but separately - in 
the space of extra dimensions. This can be described by the standard Friedmann-Robertson- Walker (FRW) metric if 
we allow for different scale factors in 3D and the UED: 

ds 2 = -dt 2 + a 2 {t) ll] dx i dx j + b 2 {t)% q dy v dy q , (7) 

where 7^ and "f pq are maximally symmetric metrics in three and n dimensions, respectively. Spatial curvature is 
thus parametrized in the usual way by k a = —1,0, 1 in ordinary space and kf, = —1,0, 1 in the UED. Of course, one 
could imagine a model that is not described by the metric J7J. For instance, there is no theoretical or observational 
argument against having separate scale factors for each extra dimension. We choose this metric because it is the 
simplest realistic alternative for studying dynamical extra dimensions. 
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With our choice of metric, the energy-momentum tensor must take the following form: 

-p 

T\ ; : ! f jPa 

7> 6 



which describes a homogeneous but in general anisotropic perfect fluid in its rest frame. The pressure in ordinary 
space (UED) is related to the energy density by an equation of state p a — w a p (jpb = Wbp)- 
The nonzero components of the field equations l@J are then given by 



3 I - I +3-^ + 3n-- 
a z a b 



a b n(n — 1) 



2-+ - 

a \ a 



k„ b ab nin — 1) 
-f +n- + 2n--+ 1 ' 
a 2 - b ab 2 



3- + 3 - + 3^ + (n-l)~+3(n-l)~ 



a 6 (n-l)(n-2) 



6 2 
6 2 

h 
b 2 



K 2 p, 



= -K 2 W a p, 



-K 2 WbP, 



(9a) 
(9b) 
(9c) 



where an overdot denotes differentiation with respect to cosmic time t. From conservation of energy T A . A — we 
find, furthermore, 



P 
P 



-3(1 + W a ) n(l + W b )-r . 

a b 



For constant equations of state this can be integrated to give 

-3(l+w <> 



P = Pi\ — 



-n(l+w b ) 



(10) 



(11) 



We will use a subscript i to indicate arbitrary initial values throughout. 



B. On energy density and pressure 

The energy density and pressure appearing in the above equations are not the usual three-dimensional quantities but 
their higher-dimensional analogues. The pressure in some direction X A is conventionally defined as the momentum 
flux through hypcrsurfaccs of constant X A . This can be expressed as 



PA , E 



If 2 V 

K A \ _ 



■|^/(k,x,i)d 3+ "k, (12) 



where k^ is the momentum in direction X A , g is the statistical weight and /(k, x, t) gives the phase space probability 
distribution. Isotropy in our model means that p a = (k 2 /i?) = (kf/I?) = (k§/£j) and pi, = (k|/_E) = ... = 
(kl +3 /E). Therefore we find 



3p a + np b = p- (— J (13) 

where p = (E) and m is the mass of the particles producing the pressure. In the case of different particle species one 
has to sum over all of them in Eq. 112fl . and the mass appearing in Eq. (|13fl can then be interpreted as the effective 
mass of all particles. For highly relativistic particles, Eq. p3|l reduces to 

3w a + nwb = 1. (14) 

Setting w a — Wb (corresponding to a completely isotropic (3+n+l)- dimensional universe) would then result in the 
equation of state 

V=^—- (15) 
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As expected, for n — we find the well-known relation for a relativistic gas in (3+1) dimensions. 

How can we recover standard cosmology with this setup? Let us first consider the case of highly relativistic particles 
and static, compact extra dimensions. Equations (|9a(l and (|9b|l are then equivalent to the ordinary Friedmann 
equations with three-dimensional energy density p^ , pressure , and an effective cosmological constant A e ff given 
by 



P 



(3) = 
(3) = 



P 
A e ff 



n(n — 1) fa 
2 ¥ 



(3) 



,-3(l+w ) 



(16a) 
(16b) 

(16c) 

For vanishing extra- 



Moreover, from Eq. Ijlll) we then find the standard cosmological evolution of p 
dimensional curvature, all we need in order to recover the familiar case of (3+1) -dimensional radiation domination 
is to set w a = 1/3. However, this forces us to allow for different pressures in ordinary space and the UED, since 
according to Eq. (|14|) Wb must then be close to zero. That is, the extra dimensional pressure (and momentum) must 
be negligible, and thus the KK tower must be largely unpopulated compared to the 3D radiation in order to reproduce 
standard cosmology for a radiation-dominated stage with w a ~ ^ ■ 

We also note that if LKPs do form a substantial part of the dark matter, then from a (3 + n+l)- dimensional point 
of view (ignoring a possible epoch of vacuum energy domination) the universe is always dominated by relativistic 
particles. This is because any standard model particle with extra-dimensional momentum is automatically relativistic 
since m -C l/R ~ 1 TcV. During ordinary radiation domination, on the other hand, the contribution from dark 
matter is negligible - but then the dominant part of the energy density is relativistic anyway. Therefore, Eq. (|14H is 
always valid in our model. Of course, we still want w a ~ in order to describe what looks like a 3D matter dominated 
universe, so we have to set 



1 

w b - 

n 



(17) 



in that case. 

On the other hand, for a time-dependent scale factor b, Eqs. I|9a|) and (|9b|l can still be cast in the standard 
cosmological form by absorbing all terms containing factors of b and its derivatives into an effective three-dimensional 
energy density and pressure: 1 



PcS 



„(3) 
PcS 



M 2 Pl 



K 2 p 



2>n 



a b n(n — 1) 



(3) 



8tt 



k w a p + n- + 2n— - + 



h 



a b n(n — 1) 



fa 
IP 



(18a) 



(18b) 



Note, however, that the actual three-dimensional energy density does not evolve in the standard manner: 

p( 3 ) = V (n) p cx a -^ +w -h- nWb , (19) 



(3) 

and there is no reason to expect that p^ ff would either, if at the same time we want to keep the standard behavior 
of a. Finally, an era of effective radiation (matter) domination corresponding to a ~ t 1 !" 1 [a ~ t 2 / 3 ) and w e g ~ 1/3 
(w e ff ~ 0) need not correspond to actual radiation (matter) domination, i.e. w a 3> Wb (wb 3> w a ). 

Finally, we would like to mention that according to Eq. © one expects the gravitational coupling constant Mp l 
to vary with a time-dependent V/ n ) ~ b n . Since in the UED model all particles are allowed to propagate in all 
dimensions, a similar case can be made for other interactions. Therefore, any nonstatic solution for b must obey the 
tight observational bounds on the allowed cosmological variation of the gravitational and electromagnetic coupling 
constants (see, for example, [35|, |3f| |3?J and references therein). 



1 For an alternative definition of the effective pressure, see l34l . 
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III. SOLUTIONS WITH STATIC EXTRA DIMENSIONS 




Static extra dimensions is the only case considered so far in the UED context and we saw above that there may be 
severe problems reproducing standard cosmology otherwise. Let us therefore study whether the field equations © 
admit static solutions for 6. Taking the difference of twice Eq. lj§c)l and the sum of Eqs. (|9af) and three times (|9bl) 
gives 

5 + 3jJ + (n-l)^) +(n-l)£ + "-,;■;■ ^ p = 0. (20) 

From this we can immediately read off a necessary condition for exactly static extra dimensions: 

(n + 2)(n-l)^ = (l-3w a + 2w b )K 2 p. (21) 

If the extra dimensions are flat (for n = 1, the curvature is automatically zero), this requires the universe to be 
empty (p — 0) or the equations of state to satisfy the following constraint: 

1 - 3w a + 2w b = , (22) 



which agrees with the five-dimensional models considered in 38, 39] . In both cases, setting 6 = 6 = reduces Eqs. I|9a|) 
and (|9b|) to the ordinary Friedmann equations for a (3 + 1) - dimensional universe and Eq. I|9c|l to a linear combination 
of these. The static solutions are therefore consistent with the full set of field equations. The particular combination 
w a = 1/3, Wb — has been found and verified before [23|, |2(| . For flat extra dimensions, there are thus two ways 
of getting static solutions - although the case of an empty universe is, of course, not particularly interesting. 

If, on the other hand, the extra dimensions are curved, Eq. I|21|l requires p to be constant for static 6. 2 Unless a is 
also static, Eq. llOf) then implies w a = —1- The origin of such an energy density could, for example, be a (3+n+l)- 
dimensional cosmological constant A, for which p = A/n 2 and w a = Wb = — 1. Setting 6 = 6 = then reduces Eqs. © 
to the ordinary Friedmann equations for a (3 + 1) -dimensional de Sitter universe, with an effective energy density 

„(3) _ M pi ( a n(n-l) k b \ _ M 2 pl A 

PcS ~ ~~&T \ 2 ¥J ~ 47r(n + 2)' {23) 

where the second equality follows from Eq. 1|21|) . Curved, static extra dimensions are thus possible in principle, but 
only for constant p. 



IV. RADIATION DOMINATION 



Recent measurements of the cosmic microwave background indicate that the universe is flat to a high degree of 
accuracy We therefore set the 3D curvature to zero. Moreover, we have shown that static extra dimensions are 
incompatible with extra-dimensional curvature in the case of non-negative pressure in 3D. Although the latest results 
from type la supernovae observations strongly suggest that the presently dominating energy component indeed does 
have negative pressure |4lj . such a component, be it a cosmological constant or a quintessence field, is believed to 
be negligible until relatively recently and it therefore cannot provide static, curved extra dimensions at earlier times. 
Thus we take k a — kb = from here on. 

Shortly after freeze-out the LKP contribution should be negligible and the total energy density of the universe is 
therefore dominated by ordinary radiation, i.e., the KK tower is largely unpopulated compared to relativistic particles 
with no extra-dimensional momentum. Equation (|14fl then implies that the extra-dimensional pressure should be 
negligible. Now, combining Eqs. (|14fl and (|22|l we find that w a = 1/3, u>& = is the only choice for the equations of 
state in a universe that is dominated by relativistic particles and has exactly static UED: 

b(t) = h . (24) 



2 Strictly speaking, this is true only for constant w a , which implies a constant iuj according to Eq. 1141 . From standard cosmology, 
however, we expect long periods of approximately constant w a (for example, during matter domination). 
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Furthermore, we have already seen that this choice reproduces standard cosmological radiation domination in that 
the scale factor and energy density evolve as 

a(t) = a, (j) 2 , (25) 

p(t)= Pl {^y\ (26) 

Thus in the UED scenario, static extra dimensions arise naturally during radiation domination. 

However, from our discussion on pressure, we do not expect w a and w b to take exactly these values and we must 
therefore consider perturbations of the static solution. So let us assume that ^ wj, < 1 and look for solutions of the 
form 

a(t) = a*(t) + 8 a (t) , (27a) 
b(t) = K+S b (t), (27b) 
P (t) = p*(t)+S p (t), (27c) 

where a star denotes the unperturbed solutions (|2~H) - l(2!))l and |(5 a |/a*, \8 b \/b*, \S p \/p* -C 1, with corresponding relations 
for the time derivatives of these quantities. The linearized versions of the field equations @ are then given by 

— 1 + 2- 2— \+n— — = —p, (28a) 

a* / \ a* a* / a* o* 3 



2^1 


; 1+ 


da 




a» 1 




a* 


a* 


_Sa _ 




) + 


(a* 


a* 


a* 







5a 




a* 













(l + 2 ^_ 2 M +n ^ + 2n ^ = --(l-n Wb )p, (28b) 
a* / \ a* a* y a* 3 

a* / <5 n <5„ \ / a* \ z / 5„ 8 n \ n — 18 b , , d» 8 b k 2 , 

a* \ a* a* y \ a */ V a * a * I 3 o* a* o* 3 

Subtracting the unperturbed equations this can be rewritten as 

a*\ 2 / S a S a \ d» 5b _ k 2 

6* 3 



a* / \ a* a* / a 



a* \ a* a* y \ a */ V a * a * J & 

— +3 — — = k w b p* . (29c) 

Equations l|29b|) and i|29cll are uncoupled differential equations for S a and S b which can be solved after inserting the 
unperturbed solutions (|25(l . (|26|l . Equation (|29a|l can be interpreted as defining 8 p and is thus automatically satisfied. 
The general solutions are found to be: 



5 a (t) = Ail-J + A 2 [-\ - -nK 2 Pl t 2 ai [-) J Wbix)^- 1 - r 1 ) dx , (30a) 

/.x-l/2 -t 

8 b (t) = B 1 +B 2 [--\ +2n 2 pit 2 b t J Wbix^x' 1 -t- 1/2 x- 1/2 )dx, (30b) 



where A\Ai-,B\,Bi are integration constants fixed by the initial conditions at time i,-. Note that and B\ can just 
as well be regarded as part of the initial conditions for the unperturbed solutions a» and &* respectively. Thus for 
w b = we find only decaying solutions for the perturbations in both scale factors - i.e. the static solution given by 
Eqs. (|24|l . (|25() is stable under small perturbations, as claimed previously (without proof) in |40l |. 
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In the case of Wb 7^ 0, however, there might also exist a growing solution. For example, a constant Wb gives 3 

s a(t) ~ -^w b riK 2 pitfai r^-J In f ^ , (31a) 
~ 2w b n 2 Pl t 2 l b l In . (31b) 
Since the relative perturbations of the scale factors grow only logarithmically in both cases, 

i^l, ^ acini, (32) 
these perturbations can still be expected to remain relatively small during radiation domination after LKP freeze-out. 

V. MATTER DOMINATION 

We have shown that there are approximately static solutions for the UED during radiation domination. In fact, 
we found that w a — 1/3 and w b = is the only possible choice for the equations of state that can give exactly static, 
fiat extra dimensions in a universe dominated by relativistic particles. However, as we remarked before, the universe 
is always dominated by relativistic particles if a significant amount of the dark matter is made up of LKPs. It is 
therefore clear that there are no exactly static solutions for the UED during matter domination, i.e. w a ~ 0. 

Having excluded exactly static extra dimensions, the next case of interest would be slowly evolving solutions: 

(33) 

But what ansatz should we make for a? A significant, long-term deviation from the usual time evolution during 
matter domination would most likely alter the predictions of standard cosmology concerning, e.g., large-scale structure 
formation. Let us therefore first examine whether there are any solutions to the field equations that give 

2 

a(t) = o*(t) = ai (£\ 3 . (34) 

We noted before that w a = implies w b — 1/n. The field equations I© can then be rewritten as two homogeneous 
equations for a and b and one defining equation for p. With the above expression for a we get 

a *\ 2 , a+b n(n - ,. )r 

n ^& + — e - UJ = T p > (35a) 

, (35b) 

^- + 7 = 0. (35c) 
in a* 

Differentiating Eq. I|35c|l . one finds that the last two equations are inconsistent with each other for n 7^ — 3, so there 
are no solutions with a(t) — a*(t) and w a — 0. Neither can there be any solutions with w a ss and a of the form 

o(t) = a*(t) + 5 a (t) , (36) 

where |<5 a |/a*, |5 a |/&*, |5 a /ffl*| <C 1. This is because, to zeroth order, inserting such an ansatz would leave the unper- 
turbed equations unchanged. Of course, if we allow for rapidly oscillating a(t) then S a and 8 a need not be small. We 
do not consider such behavior here. 



b 




a 




< 




b 




a 




A more realistic time dependence of Wi is considered in Sec. I VII where we find a more rapid growth of 8f,. 
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So are there any solutions to the field equations at all during matter domination? With the quite general ansatz 



a(t) = a>i[ — 



(37) 



we find that all solutions are of the form 



b(t) = bi 



t + B 
U + B 



(38) 



Some of these, namely 



3±y / 3n(n+2) 
3( n+3) ' 

n(n+3) • 



(39) 



are vacuum solutions. In fact, they are known as Kasner-type solutions and have been found before under the 
assumption of an empty universe and a power-law behaviour of both scale factors |42j . The only additional solutions 
appear when n ^ 1 and are given by 



2 

n+l ' 



(40) 



with B being an arbitrary integration constant. Although they describe a nonempty universe they have a static scale 
factor a. Of course, for n = 3 we knew of this solution beforehand, since w a = and Wb = h is the static, radiation 
dominated solution with a and b interchanged. 

With no suitable solutions describing actual matter domination (w a Ctut), we turn to the possibility of an era of 
effective matter domination, i.e., p^j = 0. Inserting the defining Eqs. (|18|l into the field equations @ then yields the 
familiar result a oc i 2 / 3 by construction. Using Eq. i|14|) we can write the remaining equations as 



k 2 p 



3nt 2 




4(l+ni£) 



8 + Ylnt\ + 3n(n - l)t 2 

4 2nb n(n - 1) f b\ 
3t 2 ~ + Tb + 2 I b I 



(41a) 



(41b) 



(41c) 



Equation 141afl has the general solution 



b = B 1 




l/(n+l) 



(42) 



where B\ and Bi are integration constants to be fixed by the initial conditions. This corresponds to decaying, 
bouncing extra dimensions. Although such a behavior of b might be possible, the corresponding evolution of w a is not 

(3) 

- in fact it is singular. Indeed, in our model there is no physical motivation for why jr ff should vanish for nonstatic 
b, and it is therefore no surprise that we get unphysical solutions from imposing it. 



VI. TRANSITION PERIOD 



So far we have focused on the two extreme cases of having zero pressure in either the UED (radiation domination) 
or 3D (matter domination). In a final attempt to find solutions which resemble the standard cosmological evolution 
of a (in particular during matter domination), we will now make a more general numerical study of the transition 
from an era of radiation domination with approximately static UED to one with a sizable energy density contribution 
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from LKPs. In order to do this we make the approximation that the LKPs have only extra-dimensional momentum, 
which should be valid for temperatures below ~ 1 TeV. This allows us to split the energy density and pressure into 
two parts: 



P 
Pa 

Pb 



Pr + Pra , 



w a Pm = y 



W b p r + W b p m 



(43) 
(44) 

(45) 



where r and m denote ordinary particles (radiation) and LKPs (matter), respectively. Neglecting interactions, the 
energy-momentum is separately conserved and Eq. Ijl0|l gives: 



Pr 



Pr 



Pm Pm 

Now introduce the dimensionless variable t' = t/t{ and rescale 

b 



-(n+l) 



a 
a, 



Using a prime to denote differentiation with respect to t' , Eqs. (|9 



2— + - 



b" a' . 
n— + 2n — - 
b a o 



! cj) become 
n(n — 1) 



3— +3 

a 



+ (n-1)— + 



b" 

D— + 3(n 
o 



a'b' (n-l)(n-2) 



b'\ 2 


K 2 tfp ri 


^) 


3a 4 6" ' 


b'\ 2 




~b) 


£ na 3 b n+1 



where e = p mi / p Ti , and from Eq. (|9a|l we get 

3 



2.2 

K t i p Ti 



1 



a' b' 
n—- 
a o 



n(n — 1) 
6 



(46a) 
(46b) 

(47) 

(48a) 
(48b) 

(49) 



Starting from the solutions l|27|l . (|30| l for radiation domination and approximately static extra dimensions, the 
appropriate initial conditions are given by 



o(l) 
6(1) 



, 1 n 
fl «=2-3 £ ' 



e< 1 



(50) 



Here, we keep track of terms linear in e in order to be consistent with the expected behavior of slowly growing extra 
dimensions (as opposed to the case of exactly static extra dimensions that would result from e = 0). Of course, one 
could in principle imagine different initial conditions, but the ones chosen above correspond naturally to the setup 
presented here and in Sec. IIVI 

The numerical solutions to Eqs. (|48|l are plotted in Figs. Q and |21 for different numbers of extra dimensions and 
values of e. In the beginning, we find the behavior expected from our discussion in Sec. IIVI 4 - very slowly growing 
extra dimensions and an expansion of 3D that corresponds to the usual radiation domination a oc t' 1 / 2 . However, as 
soon as the LKPs make up roughly 10% of the total energy density of the universe, the extra dimensions start to 
expand at a rate comparable to a. Such a rapid expansion of the UED is already ruled out by present bounds on the 



With the explicit time dependence of Wf, given by Eq. EU, we find that <5t, grows as t ' , which is much faster than logarithmically. 
The radiation dominated era is therefore shorter than expected from Eqs. I3UI . although the qualitative behavior remains the same. 
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log? 



FIG. 1: The evolution of the scale factors a (dashed) and b (dash-dotted), as well as the LKP energy density pm/p (solid) for n =1 
(thin), 2 (medium), and 7 (thick) with e = p mi / Pri = 10 -7 . For p m /p ^ 0.1 the extra dimensions are nearly static and the evolution of a 
reproduces the radiation dominated regime of standard cosmology to a very good approximation. For pm/p > 0.1 however, neither a nor 
b show the desired behaviour. 

time variation of the electromagnetic and gravitational coupling constants (see, e.g., [3^, |3(| |3tJ ) ■ Moreover, when 
the extra dimensions are increasing, our scale factor increases less rapidly than i' 1 / 2 - instead of approaching i' 2 / 3 as 
predicted by standard cosmology. 

The behavior of the scale factors as described above is stable against perturbations in the initial conditions of a' 
and b' of the order of e. Allowing for even larger perturbations, the only qualitatively different behavior we find is 
collapsing b and collapsing or inflating a - these solutions are obviously not viable alternatives either. 

VII. DISCUSSION AND CONCLUSIONS 

The identity of dark matter is one of the most challenging puzzles in modern particle physics and cosmology. 
Recently, it has been noted that models with universal extra dimensions provide a natural WIMP candidate that 
could make up a significant amount of the dark matter today. This subsequently led to a great deal of interest 
in studying the detection properties of these particles. However, the estimates for today's LKP abundance depend 
crucially on the underlying assumption of static extra dimensions and a standard cosmological evolution history of 
the universe since the time of freeze-out of these particles. 

In this article we have studied in detail whether one can expect such behavior without adding an explicit stabilization 
mechanism. To this end we have analyzed cosmological solutions to Einstein's field equations in (3 + 71+1) dimensions 
that are appropriate to describe a universe with UED. More specifically, our setup is given by a homogeneous FRW 
metric with two different scale factors and the assumption that LPKs make up the dominant part of the dark matter, 
i.e., the universe is always dominated by relativistic particles. 

We find that a natural - and in fact the only - way to get exactly static extra dimensions in this scenario is to 
set w a — | and wi, = 0. This also reproduces the usual radiation-dominated behavior of the scale factor in 3D. 
Allowing for + 1 Wf, <C 1, which is much more realistic in the UED scenario, we still find approximately constant 
extra dimensions and a ~ i 1 / 2 . 

However, during matter domination (w a — 0) there are no static solutions for the extra dimensions. Even worse, 
there are no solutions at all that are consistent with the standard matter dominated behavior of the scale factor in 
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FIG. 2: This plot shows the evolution of the same quantities as in Fig. 1, this time with e = 10 3 . As expected, a larger value of e gives 
a shorter period of radiation domination. 

3D, a(t) oc t 2 / 3 . With a more general ansatz a(t) oc t x we do find solutions, but for x ^ they all describe an 
empty universe. Demanding u> c ff = instead of w a — we get the usual behavior of a{t) by construction, but the 
corresponding solutions for b(t) and p(t) are unphysical. 

In the reasonable approximation that the LKPs have only extra-dimensional momentum, we have also performed a 
numerical analysis of the transition from a radiation dominated universe with approximately static UED to one with 
a sizable energy density contribution from LKPs. The evolution of b is generically found to be much too rapid given 
the present bounds on the time variation of the electromagnetic and gravitational coupling constants, and a does not 
show the standard behavior either. 

To summarize, we have shown that, within our framework, an explicit mechanism is needed not only to stabilize the 
UED but also to reproduce standard cosmology in 3D during matter domination. Although one could consider more 
complex models, e.g., with different scale factors for each extra dimension, we believe that finding static solutions - 
or indeed any solutions which give standard cosmology for both radiation and matter domination without obviously 
violating experimental bounds on the evolution of the extra dimensions - is a generic difficulty of this scenario. 
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